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Abstract 

We discuss asymptotics of the number of states of Boson gas whose Hamiltonian 
is given by a positive elUptic pseudo-differential operator of order one on a compact 
manifold. We obtain an asymptotic formula for the average of the number of states. 
Furthermore, when the operator has integer eigenvalues and the periodic orbits of pe- 
riod less than 2tt of the classical mechanics form clean submanifolds of lower dimensions, 
we give an asymptotic formula for the number of states itself. This is regarded as an 
analogue of the Meinardus theorem on asymptotics of the number of partitions of a pos- 
itive integer. We use the Meinardus saddle point method of obtaining the asymptotics 
of the number of partitions, combined with a theorem due to Duistermaat-Guillemin 
and other authors on the singularities of the trace of the wave operators. 



1 Introduction 

The subject on asymptotics of the number of partitions p{E) of a positive integer E is one 
of the main topics in classical analytic number theory, and there are many results in this 
area since the celebrated work of Hardy- Ramanuj an ( [HR2] ) appeared. As a simplest form 
of their theorem, one has 

p(E) ^eH-J , E ^ oo. (1) 

aeVs 

There are many other works on the asymptotics of p{E) or its variants. For example, after 
Hardy- Ramanuj an's work, Rademacher ([Raj) gives an expansion of p{E) in a convergent 
series. See |An] . [EL], |Rilj . |Ri2j for some of them. 
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The problem which we are going to address here is on asymptotic behavior of the number 
of states of Boson gas which is, as in the title, regarded as a spectral analogue of the number 
of partitions. 

In this section, we explain the problem which we consider in this paper, and we state the 
main results after we give some accounts on the motivation for the problem. The method we 
use to obtain asymptotics of the number of states of Boson gas is Meinardus' suddle point 
method ([Xnj, [Mi]) for obt aining asymptotic formula of the number of partitions. This 
method has some applications to mathematical physics, such as finding asymptotic growth 
of the state density for p-branes ( [Ac] . |BKZj . |EORBZj ). We review, also in this section, 
Meinardus' theorem and explain roughly the method. 

1.1 A problem in quantum statistics 

Let P be a first order strictly positive elliptic pseudo-differential operator of order one on 
a compact connected smooth manifold M of dimension n with principal symbol < p G 
C°°{T*M \ 0). For each positive integer N, consider the operator 

^ (i) 
= ^ / ® ■ ■ ■ ® / ® P ®/ ® ■ ■ ■ ® / (2) 

i=l 

on the L^-space L'^{M^) on the iV-fold product space = M x ■ ■ ■ x M (iV-times). We 
have the natural action of the symmetric goup &n on L'^{M^), and this action commutes 
with the operator . Therefore, the operator can be restricted to the symmetric and 
anti-symmetric subspaces 7-^5^, (see Section [2] for the precise definition). The restrictions 
P^ = P^ly^N and PP = P^l-j^N are regarded as the Hamiltonian for the system of N- 
particles consisting of Boson and Fermion, respectively, each particle of which is dominated 
by the Hamiltonian P, and it is assumed that there are no interaction among each particle. 
In this paper, we restrict our attention to the systems of Boson. Let Spec(P^) denote 
the spectrum of the Bose Hamiltonian, P^ . It is easy to see that Spec(P^) = Spec(P^), 
and hence Spec(P^) consists of eigenvalues of finite multiplicities. But the multiplicities of 
each eigenvalues of the Bose Hamiltonian P^ are different from that of P^ . Denote the 
multiplicity of the eigenvalue E G Spec(P^) by Q{N, E). The quantity VL{N, E) is naturally 
regarded as a number of states of A^-Boson system in the energy E. 

Natural question here is to find asymptotic formulas for VL{N, E) when E and get large 
and E and A^ are related by a relation, like E/N const. But, it seems to be difficult to 
find them. Thus, it would be also natural to consider the asymptotics as A^ — 00 of the 
density of states, which in this case would be given by 

^ Vl{N, E')5e' or more specifically ^ n{N, E') 

E'&Spcc{P^) E'e[E,E+eN] 

with a suitable sequence ejy > 0, however it is not so transparent how to choose e^. These 
two problems would be rather natural in physical points of view, because it is important to 
consider systems consisting of large but finite number of particles. Although there are some 
works (see jEcj ) in this direction, it seems to be hard to resolve. 
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Thus, the next natural quantity to consider is to take an "average" of the quantity 
Q{N, E) in the number of particles, A^. So, we define 

n{E) = J2 ^(^' E), Ee r(p), ^](o) = i, (3) 

N>1 

where we set 

r(P) = {0} U U Spec(P^) C M>o. (4) 

Af>l 

It is easy to see that the set r(P) is discrete subset of ]R>o without finite accumulation points 
and the sum in ([3]) is finite for each fixed E G r(P). We call the quantity VL{E) the number 
of states of Boson gas. Then our problem is to find asymptotic formulas for VL{E) as energy 
E getting large. 

Before stating our result on the asymptotics of f^(-E), we should mention why this kind of 
problem can be regarded as an analogue of the asymptotics of the number of partitions p{E) 
for positive integer E. To explain it, consider the generating function Gp of the quantity 

Gp{r)= J2 my-""^ (5) 

E&T{P) 

which is called the grand partition function for Boson gas. Note that the usual grand partition 
function has one more parameter, called chemical potential. The chemical potential plays 
an important and essential role in the usual setting-up for analyzing, for example, the Bose- 
Einstein condensation. The main reason for removing the chemical potential is that we 
assumed here that the operator P is strictly positive, and so a phenomenon like Bose- Einstein 
condensation does not occur. 

Now it is easy to explain the relation between the number of states Vt{E) of Boson gas 
defined for E G r(P) and the number of partitions p{E) for positive integers E. Namely, we 
have the following: 

oo 

Gp{r) = - e-'n-\ (6) 
1=1 

where < Ai < A2 < ■ ■ ■ denotes the eigenvalues of P counted repeatedly according to 
their multiplicities (see Lemma [2.3p . Thus, if = £ for each £ and the multiplicities of the 
eigenvalues are all one, then Vt[E) is nothing but the number of partitions. 

We mention here about the motivation to address the problem stated in the above. One 
of main motivations comes from the paper [TZj on the various aspects of asymptotics of the 
multiplicities m7v(A,/i) of the weights /i in the tensor power V®^ of a fixed irreducible 
representation V\ with the dominant weight A of a compact Lie group. By the Borel- 
Weil theorem, the irreducible representation V\ can be realized as the space H^{0\, Lx) 
of holomorphic sections of a suitable line bundle Lx over the coadjoint orbit Ox through 
A. Then, as a representation, V®^ is naturally identified with H^(Ox,L^^), where Ox is 
the A^-times product space of the coadjoint orbit Ox and is the A^-th external tensor 
power of the line bundle Lx- In this point of view, it would be natural to think Ox as a 
classical phase space of a particle, and hence Ox would be regarded as a phase space of a 
system of A^-particles. Then, the multiplicities m^{X, fi) would be regarded as the number of 
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stationary states for the action of the maximal torus. Therefore, the problem on the number 
of states Q{E) of Boson gas would be regarded as an analogue of the multiplicities mjv(A, /i) 
in this context. (More precisely, one should consider the multiplicities in the symmetric 
tensor power Sym^(VA) for a complete analogy of 'Boson' gas.) 

1.2 Results 

First of all, we state a theorem which is valid without any assumption on the eigenvalues of 
P. In such a general case, it seems to be still hard to obtain the asymptotics of the quantity 
fl{E). Thus, we take its average. Define the quantity D{E) by 

D{E):= Yl ^(^)- (7) 

Ler{P)n[o,E] 

Theorem 1 We have 

D(E) = e^"^"+' , r(P)9E^oo, (8) 
where the constant Bn is given by 

In particular, we have 

lim E'^ log D( E) = Bn. (10) 

r(P)3£;^oo 

In Theorem [H we do not need any assumption on the spectrum of P nor classical me- 
chanics. One may think Theorem [1] as an analogue of the usual Weyl asymptotic formula 
for the operator P itself. 

To state our next results, let us prepare some notation. Let 

^p(s) = ^A7^ sec (11) 

denote the spectral zeta function. It is well-known ( |DGj . |Shj ) that the series Zp{s) con- 
verges absolutely on the right-half plane Re (s) > n, and it is meromorphically continued to 
the whole complex plane. The poles of the function Zp{s) are located at 

s = n,n-l,...,l,-l,..., (12) 

and it is holomorphic at s = 0. We set 

Kj = Aj({n - j + l)T{n - j), j = 0, 1, . . . , n - 1, Aj = ReSs=n-jZp{s), (13) 

where ( and F denote the Riemann zeta function and the Gamma function, respectively. 
For example, we have Aq = (27r)~'^Vol(S) with S = p~^{l) C T*M \ being a level surface 
of the principal symbol p and Vol(S) being the volume of S with respect to the Liouville 
measure. Thus, we always have Kq > 0, but, in general, Kj could be non-positive (but it is 
real). Let : S ^ S denote the Hamilton flow generated by p (restricted to S). 
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Theorem 2 Assume that the eigenvalues of P are integers. Assume also that, for each 
period T with < T < n of a periodic orbit of the classical Hamilton flow (ft, the fixed 
point set of ipx is a union of finite number of connected submanifolds and each submanifold 
is clean and of dimension strictly less than dimS = 2n — 1. Then, we have the following 
asymptotic formula: 

n{E) = CE^ explx^'^fixEm + 0{E-^')), (14) 

where the positive constant xe depending on the energy E satisfies 

1 

"^\l + 0{E-^)) (15) 

{precise definition of the positive number xe is given in Section^, and the function f{x) on 
the positive real numbers is a polynomial given by 

n-l 

/(x) = J](n-j + l)i^,x^-. (16) 
The constants C, k and Ki in ([H]) is given by 

l-2Zp(0) 

C = jnKo) , 

det(P)i/27r(n+l) ^ ' 
_ Zp(0)-l-n/2 

where det(P) is the zeta regularized determinant of P and the constant ki is any number 
satisfying < Ki < when n > 2, and < Ki < 1/4 when n = 1. 

Recall that a submanifold Z in E of fixed points of a diffeomorphism : S — S is said to 
be clean ( |DGj ) if, for each z & Z, the set of fixed points of the differential d(pz '■ T^Tj T^S 
equals the tangent space of Z at z. A reason why we assume a cleanness condition for the 
classical Hamilton flow in Theorem [2] is that we need to use a theorem in [DGj on singularities 
of the trace of the wave operator e~**^. See Sections II ■4[ 13.21 

In Theorem [2l the assumption that the eigenvalues of P are integers implies that the 
classical Hamilton flow Lpt on T*M is periodic with common period 2tt. This follows from 
Egorov's theorem. But, there might be other period less than 27i. The assumption that 
the flxed point set of (pT for a period T of a periodic orbit is of dimension strictly less than 
2n — 1 implies that 27r is the least common period of the flow. There is a technical reason 
why common periods less than 27i are excluded. See Section 13. 2[ 

One of examples of the operator P satisfying the assumption in Theorem [2] is the operator 
P = a/ Acan + {n — 1)^/4 + {n — l)/2 with the standard Laplacian Acan on the n-sphere. 

Note also that, according to a result of Colin de Verdiere ( [CdV] ). if the classical flow 
is periodic with the least common period 2tt, and the square P^ of the operator P is a 
differential operator, then there is a constant a and a pseudo-differential operator Q of order 
— 1 such that the eigenvalues of the operator P -\- Q — a/4 are integers. Therefore, we can 
apply Theorem [2] for these systems if the cleanness assumption is satisfled. 
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Compared with the usual Weyl asymptotic formula for one particle, the asymptotic for- 
mula for the number of states Q{E) of Boson gas is rather complicated. In particular, all the 
positive poles affects the exponent of the asymptotics whereas only the largest pole affects 
in the usual Weyl asymptotic formula. As is shown in Section SI the positive constant xe 
satisfies the estimate f[T^ . Therefore, if one takes the logarithm and then the limit to reduce 
information, one obtain the following. 

Corollary 3 Under the same assumption as in Theorem^ we have 

lim \ogn(E) = Bn, (17) 

r(P)9_B-^oo 

where the constant Bn is given in ([9]) . 

Corollary [3] is basically obtained by Knopp ( [Knj ) . where the sequence {A„} of positive 
integers is assumed to be strictly increasing, that is, all the eigenvalues of the operator P 
are of multiplicity one. Thus, one can think Corollary [3] as a tiny generalization of Knopp's 
result. One might expect that the limit formula (ITTl) might hold for general system without 
the assumption that the eigenvalues of P are integers. However, in general, since VL{E) does 
not have monotonicity in E, it might be hard to obtain it. 

In Theorem [2], the positive constant is defined as a zero of a polynomial determined by 
the residues of the positive poles of the spectral zeta function. Thus, it is not so transparent. 
However, when n = dimM = 2, it is replaced by a more concrete constant i/e as follows. 

Theorem 4 Let n = dimM = 2. Then, under the same assumption in Theorem^ we have 
the following. 



Zp(0)-2 

n{E) = CE^^ exp 



2 

l + 0(^-'^^)) 



vl ' Ve 12Ko_ 
for any < ki < 1/3, where the constants C and ue is given by 

_ / C(3)Vol(S) \ _ 1 / C(3)Vol(S) \ 
" V 27r2i^ J ' ^ " det(P)v^ V 27r2 J 

1.3 Meinardus' theorem 

Among results on asymptotics of the number of partitions, a theorem due to Meinardus 
|Me] gives a natural generalization of the above formula ([1]) of Hardy- Ramanuj an. Let us 
state a result of Meinardus here. Let {a^j^i be a sequence of non-negative numbers (not 
necessarily integers) such that the infinite product 



G(r) = n(l-e 



£>1 



converges locally uniformly for Re (r) > 0. Then the infinite product G{t) defines a holo- 
morphic function on the right-half plane. Thus we can consider the Taylor expansion of the 
function G{z) = G{t) [z = e~'^) at z = 0: 

oo 

G{z) = I + ^r{E)z'^ , z = e-\ Re (r) > 0. 

E=l 
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The Meinardus theorem is an asymptotic formula for the coefficients r{E) as E —>■ oo. 
When a„ = 1 for all n, usual number of partitions p{E) coincides with the quantity r{E). 
Associated to the sequence {a„} is the Dirichlet series 

oo 

To state the Meinardus theorem, we need the following assumptions: 

(1) There exists a positive number a such that L{s) converges for Re (s) > a. 

(2) There exists a number Cq such that < Co < 1 and the Dirichlet series L{s) is 
continued meromorphically on the set Re (s) > —Cq. 

(3) The function L{s) has just one pole at s = a in Re (s) > —Cq, and it is simple with 
residue A = ReSs=aL{s) > 0. 

(4) There exists a number R> such that |i^(s)| = 0(|lm for Re (s) > —Cq. 
(H) Define the function 9{t) (Re (r) > 0) by 

oo 

e{r) = Y,a,e-'\ (19) 

e=i 

Then, for every sufficiently small e > 0, there exists a constant C > such that for 
every t = x + iy with | arg(r)| > 7r/4, \y\ < n, we have 

Re{e{x + iy)) -e{x) < -Cx-'. (20) 

Note that, by the assumption (1), L{s) is holomorphic on Re (s) > a. The assumption (4) 
is to guarantee the convergence of an integral. Then, the Meinardus theorem is stated as 
follows. 

Theorem 5 (Meinardus ([Me])) Assume that the sequence {an}'^^^ of non-negative num- 
bers satisfies the above conditions (1) - (4) and (H). Then, we have the following asymptotic 
formula for r{E) . 



r{E) = CE'^exp 



" + ^ ^ E^{Ar{a + l)C(a + 1))^ 



X (l + 0(E-"i)), 



a 

where A > is the residue of L{s) at s = a, and the constants C , k and ni are given by 

C = -^f=(Ar{a + l)ac^ + l))W 

_ L(0)-l-a/2 

< 5 < i 

Meinardus' method of obtaining the above theorem is based on the saddle point method, 
and thus it would be natural to expect that Meinardus' method would work well for more 
general problem. Roughly speaking, Meinardus approximates the function logG(r), which 
plays a role of the phase function, by a very simple function and then estimates the error 
terms by using the assumption (H). 
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1.4 Remarks and comments 



(1) It would be clear from the statement of Theorem [2] that the spectral zeta function 
Zp{s) plays the same role as the Dirichlet series L{s) in Theorem [5l Furthermore, the 
eigenvalues of P are assumed to be integers. Therefore, the situation in Theorem [2] is, at 
least formally, contained in that of the Meinardus theorem with {a„} being the multiplicities 
of the eigenvalues. However, there are something to be careful. One of them is the fact that 
the spectral zeta function could have many positive real poles, although it is assumed, in 
Theorem O that L{s) has only one positive pole. Another is the condition (H) which is 
assumed in Theorem O In Theorem [21 there are no assumption similar to the condition (H). 
In our case, the condition (H) is a condition on the theta function ('heat trace') 

oo 

0{t) = J2 e"^'", Re (r) > 0. (21) 

1=1 

Then, the cleanness assumption for the classical mechanics in Theorem [2] makes us to be 
able to prove the condition (H) for the function 9 (see Section [3]). 

The assumptions in Theorem [2] exclude the existence of the common period less than 27r. 
This is necessary for the validity of the condition (H) for the function 6 defined above. In 
fact, if all the eigenvalues are even integers, then we have 6{x + in) = 9{x) for x > 0. This 
shows that the condition (H) does not hold for the system having eigenvalues all of them are 
even integers. By Helton- Guillemin theorem ( [Gu2j ) . if tt is a common period of the classical 
Hamilton flow, then (assuming that P being the square root of the Laplacian), the set of 
cluster points of differences of eigenvalues of P, {Aj — A^}, is just 2Z. Thus, a perturbation 
of P could have eigenvalues all of them are even integers. Thus, if vr is a common period 
of the classical flow, then, one might need to impose some additional assumption on the 
eigenvalues of P to ensure the validity of the condition (H). 

(2) As mentioned, we assumed that the operator P is positive definite, and hence a phe- 
nomenon like Bose-Einstein condensation does not occur. However, when the operator P 
has the zero eigenvalue, or when we consider a family of positive definite operators Pt with 
the least eigenvalues \i{t) which tends to zero as t — oo, then the number of particles in the 
zero energy state, or small energy states, could be large. This phenomenon is Bose-Einstein 
condensation, and in such a case, we need to use the full grand partition function defined as 

S(r,^) = n(l-e"~'")"'' 
1=1 

where the parameter /z, which is called the chemical potential, should be negative. In this 
case, the quantity Vt{E) should be replaced by 

N>1 

Thus, it might be interesting to find joint asymptotics of Q{E,fi) as E oo and fi —>■ 0. 

(3) It would be hard to describe the asymptotics of Q{E) for general system. This is because, 
in general, there are no integral representation for the quantity Q{E). Although there are 



8 



some representation for Q{E) in terms of a limit of integrals (see Lemma [2.4p . we need to 
estimate the error term. In the error term, such a quantity like the difference of eigenvalues 
in r(P) appears. In general, estimating the level spacing is quite hard. However, when the 
classical system is periodic, like a ZoUl surface, one might be able to use the results in [DGj . 
jGulj . [UZj to find a similar asymptotic formula for Q{E) in such a case. 

Acknowledgments The author would like to thank to professors T. Matsui and S. Zelditch 
for their helpful comments. The author also would like to thank to the referee who pointed 
out some mistakes in the earlier version of the paper and gave valuable suggestions to the 
author. 

2 Partition function for Boson Gas 

In this section, we give precise definition and some properties of the number of states, fi(-E'), 
of Boson gas and its partition function, which are, at least formally, quite well known in 
physics textbooks on quantum statistics ( [Fuj . |Khj ). 

2.1 Occupation number representation 

As in Introduction, let P be an elliptic, strictly positive pseudo-differential operator of order 
1 on a compact connected manifold M of dimension n, and consider, for each positive integer 
N, the operator on L'^{M^) where = M x ■ ■ ■ x M {N times) defined in ([2]). We 
define the action of the symmetric group ©iv on L'^{M^) by 

{UaF){xi, ...,Xn) = F{x^(i), . . .,X^(N)), 

ae&N, F e L\M^), (xi, . . . , x^) e M^. 

This action of ©at on L'^{M^) commutes with the operator P^ . Let Tig and Ti.^ be the 
space of symmetric and anti-symmetric functions on , respectively: 

n^ = {Fe L\M'') ■,U„F = F, ae 6^}, 
n^ = {Fe L^M"") ; U^F = sgn(a)F, a G &n}- 

The spaces Tig and Ti.^ are the state space for the system of A-Boson and A-Fermion, 
respectively. In this paper, we consider the system of Boson. We denote the restriction of 
P^ on by P^. 

Let < Ai < A2 < ■ • • t oo be the sequence of eigenvalues of P. We choose an 
orthonormal basis {<^j}j>i of L?'{M) consisting of eigenfunctions of P: Pipj = Xj^Pj, j = 
1,2,.... For any multi- index a = (ai, . . . , oat) G Z^q, we define a function $q, G C°°{M^) 
by 

Then, the function is an eigenfunction of P^ , and we have 

N 

P^$, = S$„, E = Y,K- (22) 

i=l 
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From this, the spectrum of P is easily described. But, what we need to know is the 
spectrum of the Bose Hamiltonian . Let vr^ denote the orthonormal projection onto 7i^: 

= ]^ E 

ctGSjv 

Then, clearly the state space for Boson Tig spanned by the functions vr^$a, a G Z^q. The 
symmetric group & n acts on Z^q by 

aa = (a^-i(i), . . . , a^-^N)), cr E &n, a = (ai, . . . , on) e Z^g. 



Lemma 2.1 (1) tt^^q, = 7r_^$/3 if and only if {3 lies in the orbit of Gp^ through a. 
(2) Let [a] denote the orbit of &n through a G Z^q, and set = 7r_^$Q/||7r^$Q,|| . Then, 



{^[a\}[a\&z%/&N 0''^ orthonormal basis for Ti 



N 

S ■ 



Proof. First, note that we have Ua^a = ^aa for any cr and a. Now, suppose that (3 = fia 
for some fi E &n- Then, we have 

Conversely, let vr^^a = TT^^jS. Since is an orthonormal basis for L^(M^), we 

have 

(7r^$,, = ^E^*^-' '^") = ]^ttStab(a), 

where Stab(a) denotes the stabilizer of a in ©at. In particular, TCg^a is not zero, because 
tlStab(Q;) is greater than or equal to one. Thus, by the assumption that Tig = 7r_^ $q,, we 
have 

cr 

and we have tlStab(a) > 1. Thus, for some a G ©tv in the sum of the above, ( $0-/3, $a ) is 
not zero, and hence a = a (3. This shows (1). (2) follows easily from (1). | 

Since we clearly have Pg = E^\a] with E given in (I22D, the spectrum Spec(P^) of 
P^ is descrete and coincides with Spec(P^). Thus, what is important is to understand the 
multiplicities of eigenvalues of P^ . 

For any E G Spec(P^), we set 

AT 

C^(E) = {aGZ^o;E = ^A„J. 

i=i 

We then define 

^1{N,E) := ti[C^(P)/6^]. 

Then, by Lemma EH] (2), the quantity Q{N,E) equals the multiplicity of the eigenvalue E 
of the Bose Hamiltonian P^ . 
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Lemma 2.2 For any E G Spec(P^), we have 

n{N, E) = U eZ^,;E = J2 ^kXk, N = J2^A (23) 

I k>l k>l J 

Proof. For each a G Cn{E), we define non-negative integers ^^(a) by 

iT-kia) = tl{j ; = fc}- 

Then, we have 

A'' oo oo 

j=l k=l k=l 

It is easy to see that, for G Z^q, ^^(a) = for every k if and only if /5 is in the 

©TV-orbit through a, and from this fact, we conclude the assertion. | 

The expression of Q{N, E) in (l23l) is called the occupation number representation. By 
Lemma \2.2\ the sum in the definition ([3]) of Q{E) is a finite sum for each fixed E, if the 
eigenvalues of P are positive. The quantity ^1{E), with setting ^2(0) = 1, is thus a function 
on the set T{P) defined in (jlj). For each positive number A, the set T{P) fl [0,A] is finite 
because if we write E G r(P) fl [0, A] as E = J2k ''^k^k with = J2k '^fc ( l23l) . we clearly 
have AiA^ ^ E < A. Therefore, T{P) is a discrete subset in ]R>o without finite accumulation 
points. 

2.2 Partition function 

The partition function for the number of states Q{E) of Boson gas is defined, at least 
formally, by the equation (jS]). The function Gp{t) always converges and has infinite product 
representation if P is strictly positive definite as follows. 

Lemma 2.3 The partition function Gp{t) defined in ([5]) converges absolutely and locally 
uniformly for Re (r) > 0, and on this region, Gp{t) has the infinite product representation 
©. 

Proof. Set U£{t) = (1 — e^^'^'^)^^ — 1. Then we have 

< (1 - e~^^^■)~^e-^'^ 

where we set r = x + iy. Since the heat trace 9{x) = X]£>i e~'^^^ converges locally uniformly 
for X > (see, for example, [Gi]), the sum ^ (r) converges absolutely and locally uniformly 
for Re (r) > 0. Therefore, the infinite product in ([6]) converges absolutely and defines a 
holomorphic function on the region Re (r) > 0. Now, we compute the partial product of the 
infinite product in ([6]), by using the geometric series expansion for each term, as follows. 

K 

n(l-e-^-)-'= E ^K{E)e-'^\ 
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where we set 



k=l J 



For each L > 0, there exists Kl > such that L < Xkl+i- Then, it is easy to see that 
Q{E) = VLk^{E) for each E G r(P), E < L. Therefore, for positive x, we have 



n(i - e-^n-' > E ^(^) 

e=l E(^T{P),E<L 



e 



for each L > 0. Thus the partition function Gp{t), defined in ([5]), converges absolutely and 
locally uniformly for Re (r) > 0, and satisfies 



Gpix) < - e 



£=1 

for positive number x. Conversely, it is clear that for each > and E G r(P), we have 

K 



■Q(l_e-A,x)-i < ^^^^^ 



and which, with the analytic continuation, proves the lemma. | 
If the eigenvalues of P are integers, then we have 

n{E) = ^ f e^(^+^^)Gp(x + ly) dy (24) 

for any positive number x. This formula is the starting point of our analysis. In general, 
the set r(P) is not contained in the set of integers, and hence there are no such an integral 
representation for Vt{E). However, even if T{P) ^ Z, we can approximate VL{E) by integrals 
of the form flM|) . We shall mention about this formula, although this is not necessary for 
later sections. 

Lemma 2.4 For any a; > and T > 0, we have the following. 



f^e^^''-'''^Gp{x + ty)dy 



<^£M^( niin \E-E'\] (25) 

T \E'£r(P),E'^E ' 



Proof. First of all, we note that the minimum in the right hand side in f l25p exists because 
r(P) is a discrete subset of M>o without finite accumulation points. Substituting the series 
expansion for Gp{x + iy) for the integral in the left-hand side of fl25l) . we have 



^/ e-^^--^Gp{x + ry)dy~n{E)= ^(^')^^^^^^^e(-->, 

"^""^ E'&(P),E'^E ^ ' 

and this proves fl2S|) . | 
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3 Properties of the partition function 



In this section, we review some facts about the spectral zeta functions defined in ( fTTl) . and 
then, by using these facts, we derive an approximation formula for the partition function, 
which is a tiny generalization of a lemma in |Me] . |An] . We prove the condition (H) under 
the periodicity assumtion for the classical dynamical system. 



3.1 Approximation of the partition function 

It is well-known ( |DG] . |Sh] ) that the spectral zeta function Zp{s) for the first order elliptic 
positive pseudo-differential operator P has the following properties. 

(1) The series Zp{s) converges absolutely for Re (s) > n, and it is continued meromorphi- 
cally on the whole complex plane. 

(2) All of the poles of the function Zp{s) are simple and they are given by (fT2|) . In 
particular, Zp{s) is holomorphic around the origin. 

(3) For each c G M, the spectral zeta function Zp{s) has a polynomial growth in Im (s) 
locally uniformly in Re (s) > c (except on neighborhoods of poles). 

By using these facts about the spectral zeta function Zp{s) we have the following. 

Lemma 3.1 Let Kj be the real numbers defined in (fT3|) . Then, for r = x + iy & C with 
< X, \y\ < X, the logarithm of the partition function log G'p(r) has the following expression: 

n-l 

logGp(r) = ^/^,r-('^--') - Zp(0) logr + ^^,(0) + J(r), (26) 
i=o 

where the term J{r) is given by 



1 



-Co+ico 



J(r) = — / T-^Zpis)as + l)r(.) ds (27) 

for any < Co < 1. Furthermore, for t = x + iy with < x, \y\ < x, we have \ J{t)\ = 
0{x'"'^). The branch of the logarithm in ( l26l) is described in the proof. 

Proof. Proof is the same as that in |Me] , [An] , but we review it for completeness. According 
to the infinite product representation ([6]) of the partition function Gp(r), we can define the 
branch of the logarithm log Gp(r) by 

oo 

logGp{r) = -J2^og{l-e-'n, 

i=\ 

where the logarithm in the right-hand side is its principal branch. For Re (r) > and do > 0, 
we have 

e"" = — / r-T(s) ds. 
13 



By using the above and the Taylor expansion of the logarithm, 

- iog(i - e-'n = E r"'''^ ^^"^ > ^^^^ 

fc>l 

we have 

1 rao+ioo 

- log(l - e-'n = TT- / K'-r-Xis + ms) ds 
for Re (r) > and ctq > 0. Taking aQ> n and summing over £ > 1, we have 

log Gp(r) = — / T-^Zp{s)as + l)ns)ds. 

Take < Co < 1 and set s = a + zt with —Co < cr < o"o and t G M. Then, for t = x + iy 
with < X, |y| < X and s = a + it with the above region, we have 

< |r|-"e"l*l/^ (29) 

It is well-known that ('(l + s) is of order 0(|t|'^^) with a positive constant ci. By Corollary 2.2 
in [DGj . we have |-Z'p(s)| = 0(|t|'^2) with a positive constant C2- It is also well-known that we 
have |r(s)| = 0(|t|'^^e^'^l*l/^) with a positive constant C3. Thus, we can use the residue formula 
to shift the integral contour to —Co + it, t G M. The function /(s) := Zp{s)({s + l)r(s) 
has simple poles s with —Co < Re (s) <aats = n— j with j = 0, . . . ,n — 1 and a pole 
of order two at s = 0. A direct computation shows that the residue of /(s) at s = n — j 
(j = 0, . . . ,n - 1) is KjT-'^''-^') and the residue of f{s) at s = is -Zp{0) logr + Z'p{s). 
Thus, the residue formula gives the expression fl26l) with the function J(r) given by fl27|) for 
T = X + iy with x > and |?/| < x. Note that the integral J(r) is absolutely convergent for 
this region and is of order 0(x'"°) by fl2Ul) . which completes the proof. | 



3.2 The condition (H) 

Our next purpose is to prove the condition (H) on the function 9{t) defined by f l2TI) under 
a cleanness assumption of the classical Hamilton flow ipf First, we begin with a heuristic 
argument by using the heat trace asymptotics 

9ir) ~ Cor-" (30) 

whose validity is quite well-known at least for positive r. If ( l30l) would hold for r G C, 
Re (r) > as r — >• 0, the function Re6'(r) — 6'(Re (r)) in the condition (H) would be replaced 
by more simple function co(Re (r~") — Re (r)""). Then, the condition (H) would be proved, 
at least for small \y\, by the inequality 

Re((x + zi/)-") -X-" < -ix"" (31) 

for < X < However, when Re (r) > and r — > 0, we face singularities of the 
distribution 6{it), t G M, on the imaginary axis, and hence we need to take care about such 
singularities to prove the condition (H). 
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Proposition 3.2 Assume that, for each period T with < T < n of a periodic orbit of 
the flow (ft, the fixed point set of ifj- '■ ^ H is a union of finite number of connected 
suhmanifolds and each Zj is clean and of dimension dj < 2n — 1. Then, the 

heat trace 9{t) defined by fl2T]) satisfies the condition (H), namely, for every sufficiently 
small e > 0, there exists a constant C > such that, for any t = x + iy with x > and 
X < \y\ < 71, we have 

Re{9{x + iy))-9{x) < -Cx" 



-£ 



Proof Consider the distribution fi defined by 

oo 

fi = ^e-'^^* = Tr(e~'*^) = e{zt). 
i=\ 

Then, it is well-known ( |Ch] . |DGj ) that the singular support of /i is contained in the set of 
periods of the Hamilton flow of p. Let T G [— vr, tt]. If T is not a period, then /i is smooth 
around T. Since the distribution /i is the boundary value of the function Q{x + iy) as a; \ 0, 
and since the function 6'(r) is holomorphic on Re (r) > 0, 6'(x + iy) is bounded as x \ 
when y is close to the regular point T. Now, suppose that 7^ T G [— vr, tt] is a period of the 
Hamilton flow. Then, according to Theorem 4.5 in |DGj . there exists an interval around T in 
which no other periods occur, and on such an interval, there exists functions . . . , ^^(s) 

such that 



' /'+00 
/i = ^ / e~'''(^^'ciAs)ds 



(the Fourier transform of e'^'^'aj in the distribution sense) and the functions aj{s) admits an 
asymptotic expansion of the form: 

00 

k=0 

where dj is the dimension of the submanifold Zj and it is assumed that dj < 2n — 2. The 
functions aj{s) are of order \s\~^ as s ^ —00 for any > 0. 

Then, we note that the function 6(t) with Re (r) > can be written as 



E 



e-^'e'^'aj{s) ds. 



In fact, the right hand side of the above, which we denote H{t), is holomorphic for Re (r) > 
whose boundary value as a distribution is n, and which equals that of 0{t). Thus, since the 
boundary value of the holomorphic function 6 — H is zero, we have 9 = H near T (see, for 
example, |Ho2j ). We then introduce the cut-off function x such that x = 1 on [1? 00) and 
X = on (—00, 0], and we set 

9{r)=j^[F,{T)+G,{r)l 



Fj{T)= e-^'xis)e'''aj{s)ds, G'j(r) = - xis))e'' 'aj{s) ds. 
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Clearly, |Gj(r)| is bounded when a; \ and y is close to T, where we denote t = x + iy. 
Changing the variable t = xs and using the estimate < cs^^^~^^^'^, we have 

Since we assumed that dj < 2n — 2, we have x"^'^^""^-'/^"^ < cx"^/^~". Therefore, we have 

Re {e{x + iy)) - e{x) < \e{x + iy)\- e{x) < -Cqx"" + ca;^/^-" + C 

for some constants cq, c, C > 0. Thus, the estimate in the assertion follows when y is close 
to the period T 7^ 0. 

It is rather easy to handle the case where y is close to zero. Note that, by the ellipticity 
of the operator P, the function n(x, ^) on (x, ^) G T*M \ defined by 

n(x,O=mf{t>0;(^t(x,O = (a:,O} 

is lower semi-continuous strictly positive function on S. Since S is compact, the function 11 
has a minimum to > on S. Thus, there is an interval around zero which does not contain 
any positive period. Thus, we can express 9{t) as 

0(r)=Cor-" + a(r), |a(r)| = Odrl^"'^) 
with a positive constant Cq (see the proof of Corollary 2.2' in |DGj ). Then, by fl3Tl) . we have 

Re {e{x + iy))- e{x) < -Cix-"(1 + o{x)), 
which shows the asserted estimate. | 



Proposition 3.3 Fix a real number /3 such that 1 < /3 < 1 + |. Assume that the heat trace 
6{t) defined by fl2T|) satisfies the condition (H) {which is explicitly stated in Proposition \3.2^ . 
Then, there exists constants c > 0, C > and > such that 



n-1 



Gp{x + iy) exp[— Kj 

j=0 



X 



< Ce- 



for X + iy with x > , x^ < \y\ < tc . 

Proof. As in |Me] . |An] . we devide the region x^ < |y| < vr into two parts. Note that, for 
the proof of the main theorem, we need the case that < x is sufficiently small. So, we may 
devide it into x^ < \y\ < x and x < \y\ < vr. First, we consider the region x^ < \y\ < x. By 
a direct computation using Lemma 13.11 we have 



A{x + iy) :- 



n-1 



Gp{x + iy) exp[— 



X 



3=0 



n-1 



< Cexp[^ Kj{Re (x + iy)-'^""'^^ 

j=0 



X 



} + 0(|logx| 
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where C > is a constant. We set oj = y/x. Then, |a;| < 1. It is easy to show that there 
exists constants Ci, C2 > such that 



-CiX-("-^')cj2 < Re (x + - a;"("-^') < -C2X- 



Therefore, we get 

A{x + iy) < C exp 



n-l 

-C3X-"cj2 + C4Cj2 + 0(1 logxl) 



Now, we have x'^'^^ < cj^ < 1 and < n — 2(/3 — 1). Thus we have 

A{x + iy) < Cexp[-C5CJ^x-"] < C ey.^[-c^x'''+^^'^-% 

Therefore, setting e„ = n — 2(/3 — 1), we obtain the desired inequahty. Next, consider the 
case where x <\y\< n. First, we note that Re \og{Gp{x + iyj) = log \Gp{x + iy)\. By using 
(l28l) for r = X + iy, we have 



-Re log(l - e-^^("+*^)) < e"^^" cos(A^y) - log(l - e"^^") - e"^^". 

Summing over £ > 1 in the above, we have 

\og\Gp{x + iy)\ < Ree{x + iy) + \ogGp{x) -e{x). 

Since x < \y\ < tt, we can use the condition (H), that is, for each sufficiently small e > 0, 
there exists a positive constant C such that, for any t = x + iy with x > 0, x < |?/| < vr, we 
have 

Re (^(x + iy)) - 9{x) < -Cx"". 

Hence we have 

log \Gp{x + iy)\ < logGp(x) — Cx~^. 
Therefore, applying Lemma [STT] for r = x and taking e < ra — 2(/3 — 1) small enough, we have 

A{x + iy) < exp[-Cx"^ + 0{\ logx|)], 

which shows the assertion. | 



4 Proof of the main theorem 

The purpose of this section is to give proofs of the main theorems. First of all, we will give 
a definition of the positive number xe appeared in the statement of Theorem [2j 
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4.1 The positive number xe 

Define the polynomial Pe{x) with one variable by 

n-1 

Pe{x) = j)KjX^, (32) 

j=0 

where the constants Kj are defined in (fT3i) by using the residues of poles of the spectral zeta 
function. The polynomial pe is defined by the spectral data and the fixed energy level E. 
Since Kq = AoC{n + l)T{n), Ao = (27r)-"Vol(S) is positive, we have pe{0) = -tiKq < 0. 
The coefficient of x^^^ is > 0, and hence we have Pe{x) +oo as a; ^ oo. Therefore, 
there is a solution of the equation Pe{x) = in the positive real axis. Now, we have the 
following 

Proposition 4.1 For sufficiently large E > 0, the polynomial pe{x) has just one positive 
zero. We denote the zero of pe by xe {as in Theorem^. Then, the asymptotics of xe cls 
E ^ oo is given by f|T5l) . 

Proof. First of all, let > be a positive zero of the polynomial pe- Then, we have 

n-1 

o<E<Y,in-mM—\ 

j=0 

It follows from this that yE = 0{E^^) as E ^ oo. We have 

n-1 

PsiyE) = n{n + l)Koy]^^ + ^(n - j)(n + 1 - j)Kjy^~\ 

The sum in the right hand side of the above is bounded as — > cxo, and the first term goes 
to oo. (Each Kj is a real number.) Therefore, we have PEiys) > for each zero yE of pe for 
sufficiently large E > 0. Thus, for sufficiently large E > 0, the polynomial pe has a unique 
positive zero x^; > 0. We note that 

n-1 

= nKo + ^(n - j)4 = ^^o + 0{E-^), (33) 
i=i 

which shows the asymptotics (fTSl) . | 

We notice that, without any assumption on the the eigenvalues of P and classical Hamil- 
ton flow, we can prove the upper bound for Q by using Lemma 13.11 which agree with the 
limit formula in Corollary [3l 

Proposition 4.2 We have 

n{E) < cE^e^"^^(^+«(^)), B^ = {n + 1) (^C(n + l)r(n)) (34) 
as E ^ oo with a positive constant c. 
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Proof. We set 

n— 1 
j=0 

We take arbitrary M > 0. Then, by Lemma 13. we can find positive constants Ci, C2 
(depending on M) such that 

cix-'e'^^") < Gp{x) < C2X-'*e'^(") 

for any x in the interval (0,M], where we set a = Zp{0). Thus, for any E G T{P) and 
< X < M, we have 

We take a sufficiently large positive number Eq such that xe < M for any E > Eo, and, for 
-E- > -£"0, we set a; = x^;, so that the function tp{x) + Ex takes its minimum there. Then, a 
simple computation with f|T5l) shows 

+ EXij = + 0(E-^)), 

and which completes the proof. | 
4.2 Proof of Theorem [1] 

We use Hardy-Ramanujan Tauberian theorem ( [HRlj . Theorem A) which asserts the fol- 
lowing: Suppose the sequences {i?m}m=i' Wm}m=i non-negative real numbers and the 
constants A > 0, a > satisfy the following: 

(1) Em > Em-1, Em OO] 

(2) Em/Em-1 1; 

(3) G{s) := J2(^rnG~^'^^ is convergent for s > 0; 

(4) G{s) = exp[y4s~"(log(l/s))~^(l + o(l))] as s ^ 0, where is a real number. 
Then we have 

n 

An ■.= Y,am = exp[Si?:/(^+")(logi?„)-^/(i+")(l + o(l))], 

m=l 

where B is given by 

In our situation, we denote r(P) = {0 = Eq < Ei < E2 < ■ ■ ■ } ■ Then we have Em — Em-i < 
Ai, and hence Em/Em-i — > 1 as m ^ 00. Thus the condition (1) and (2) are satisfied. The 
non-negative number am in the condition (3) corresponds to the number of states Q{Em) 
of Boson gas, and hence (3) follows from Lemma 12. 3[ The condition (4) is deduced from 
Lemma [3.11 with a = n, j3 = and A = Kq = (27r)~"Vol(S)C(n + l)r(r;,). Thus, the formula 
follows from the Hardy-Ramanujan Tauberian theorem with B = Bn- 
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4.3 Proof of Theorem [2] 

First of all, we fix < < ^ (0 < < 1/2 for = 1) and choose 0<Co<l, 1</3<1 + | 
as follows. For n = 1, we set P = ^ + j and Cq = | + |//. For n = 2, we set /5 = | + ^/i 
and Co = I + |yU. For n = 1, 2, we set 5 = | — fi. For n > 3, we take 6 satisfying < 5 < 
min{|, ^, |(| — ~ f^)} set /3 = 1 + |(1 — |) and Cq = n{fi + j)- In all cases, the 

positive numbers /?, 5 and Co satisfy < 5 < min{l/2, 4/n}, 1 + | < /3 = 1 + < 1 + 

< Co < 1 and /X = min{| - 5, ^ - f }• 

From now on, we assume that the eigenvalues of P are integers. Then, we have r(P) C Z, 
and hence we can use the integral formula fl2^ for the number of states Q{E). For every 
X > 0, we have 

Q{E)=h{E) + Ri{E), 
Ri{E) = ^ [ e'^^^+'y^Gp{x + ty)dy, 

^TT Jxl^<\y\<w 

where /3 is a fixed constant as above. We shall determine a; > by using the suddle point 
argument for the integral Ii{E). By Lemma [3.11 we have 

271" J\y\<xl^ 

Then, in the integrand of the above, we think Et + Yl^=o KjT^~"- as a phase function. The 
critical point of this phase function close to Im (r) = is given by r = x^. Thus, we choose 
X = xe in (!35|) . First, let us consider the integral h{E) in (1351) . As in [Anj , we write 
tue = Exe- Then, by Lemma [3. II and the change of variable u = y/xE, we have 

I^(E) = ^mfi+Z'p(0)+(l-Zp(0))logXfl 
271 



where we have used xe = 0{E = 0{m^^^^). By using the formula (1 + iu) 

1 — iaui — ^a{a + l)u;^ + 0(|ci;|^), and 

n-l 

("-i) 



rriE = Exe = "^{n - i)KjXE 

3=0 



we obtain 



/^(^) = ^ga;^"/(x£)+Z^(0)+(l-Zp(0))logXB ^ /2(E), 



.|<4-^ ■ (36) 



n-l 



j=0 
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where the function / is defined in ( IT6|) . By the definition of the number xe and the function 
/, one has 

n-l 



j=0 



Now, it is easy to see that lul"^ = 0{m^ " ) for every j = 0, We note that 

1 — f (/? — 1) < 0. Near u = 0, obviously we have log(l + ia;) = Odtul), and, since /3 < 1 + |, 

L('/5 I') 1 — (3 1") 

we have = 0(m^"^ ') = 0{m^ ). Thus, if we set /ii = min{-(/5 - 1) - 1, Cq/^}, 



we obtain 

h{E) = 



where we set 



A(E) = [ e-^''"^^)"' du, 

J \ui\<x'^ 

By using the definition of rin{E) and ([3SD, we have 

m-ir7„(E) = (n + 1) + 0(m^'/"). (37) 
Thus, for sufficiently large E > 0, i]n{E) > 0. Therefore, we have 

Then, we easily have 



For r]n{E), by using ( l37Il . we have 



^ ^ :(l + 0(m-^/")). 



^/Vn{E) ^/(rlTT)mE 
Therefore, since = 1 + 1(1 - I) with 0<5 < min{l/2,4/n} and 

;. = mm{--5,---}, 

we obtain 



Inserting the definition niE = Exe, and using (fT5|) . we obtain 

/ Ott / 1 \ V2(n+1) 
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n + 1 V ni^o 



By fl36|) . we conclude, for the integral h{E), that 

x~"f(xp) / 1 \ l/2(n+l) 

Inserting, in the above, the formula 

i~Zp{o) _ I nKo\ 



which is deduced from fllSp . we get the right hand side of fll4p . Next, we have to consider 
the integral Ri{E) in fl35|) . By Proposition 13.31 we have 

n-l 

\Ri{E)\ < Ce-^""^" exp[ExE + J2^j'^~e''"~'^] = Ce-^""^"e"^"^("^) 

3=0 

with positive constants C, C„, e„. Since the term 6"*""^^^" decays exponentially and since 
the term e^^ /(^b) jg same as the exponential term in the asymptotics of the integral 
Ii{E), we conclude (IT^ . This completes the proof of Theorem [21 
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